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Abstract 

An explicit, detailed evaluation of the classical continuum limit of the axial 
anomaly/index density of the overlap Dirac operator is carried out in the infinite 
volume setting, and in a certain finite volume setting where the continuum limit 
involves an infinite volume limit. Our approach is based on a novel power series 
expansion of the overlap Dirac operator. The correct continuum expression is 
reproduced when the parameter ttt-o is in the physical region < mo < 2. This 
is established for a broad range of continuum gauge fields. An analogous result 
for the fermionic topological charge, given by the index of the overlap Dirac 
operator, is then established for a class of topologically non-trivial fields in the 
aforementioned finite volume setting. Problematic issues concerning the index 
in the infinite volume setting are also discussed. 

1 Introduction 



In the last few years there have been interesting developments in lattice gauge the- 
ory concerning chirality aspects of lattice fermions and topological aspects of lattice 
gauge fields. The overlap formulation of chiral gauge theories [p], |2[, which was largely 
inspired by Kaplan's domain wall proposal ||, led to the introduction of the over- 
lap Dirac operator for massless lattice fermions @]. This operator satisfies the the 
Ginsparg-Wilson relation ||, thereby providing an explicit solution to the chirality 
problem in lattice QCD (see, e.g., [|7| for reviews). 
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A key test of any lattice formulation of a quantum field theory is whether it repro- 
duces the correct anomalies for the continuum symmetries. In particular, the axial 
anomaly in QCD explains the tt° — ► 2 / y decay amplitude || and should be repro- 
duced in the lattice formulation. Another important quantity in lattice gauge theory 
is the topological charge of lattice gauge fields. The topological charge enters into 
an explanation for baryon number violation in the standard model, and for the rela- 
tively large value of the rf mass in QCD ||. Lattice definitions of topological charge 
is therefore an interesting topic of study. A necessary condition for an acceptable 
lattice definition of topological charge is that it has the correct continuum limit. 

In this paper we study the classical continuum limit of the axial anomaly and 
fermionic topological charge in lattice gauge theory with fermion action specified by 
the overlap Dirac operator. Our approach avoids the limitations of earlier approaches 
which were only valid for small (and hence topologically trivial) background gauge 
fields. In §3 we develop new techniques (first introduced in [fTA, (v4)] and [II]) for 



evaluating the classical continuum limit of the index density/axial anomaly, which 
enable us to rigorously show that the correct continuum expression is obtained (when 
< m < 2) for a broad range of continuum gauge fields. These include topologically 
non-trivial fields and more general classes of fields which may diverge at infinity or 
have singularities. 

The situation regarding the index is more delicate. In the continuum, in infinite 
volume, the index of the Dirac operator coupled to an SU(n) gauge field A(x) = 
A't{x)T a dx IJ ' on R 4 is ill-defined in general. A related fact is that the topological 
charge, 



Q(A)= I d 4 xq A (x) , q A {x) = — \e^ P MF^{x)F pf7 {x)) , (1.1) 
'r 4 32tt z 
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is not integer valued|] or even well-defined,Q in general. However, integer values for 
Q(A) are guaranteed when A converges sufficiently quickly to a pure gauge configura- 
tion at infinity (or at the singularity in the case where A is in a singular gauge). It is in 
this case that the Dirac operator coupled to A has a well-defined index equal to Q(A) 
by the Index Theorem (since in this case A^x) corresponds to a gauge potential in an 
SU(n) bundle over the compact manifold S A = R 4 U {00} [0). Similarly, in the lat- 
tice theory the index of the overlap Dirac operator D u will be ill-defined in general in 
the infinite volume setting (i.e. infinite lattice on the whole of R 4 ). One expects that, 
in analogy with the continuum case, index D u will be a well-defined finite number 
when suitable conditions are imposed on the lattice gauge field U, e.g. a requirement 
that it converge sufficiently quickly to a pure gauge configuration at infinity. A first 
step towards clarifying this issue is carried out in §4, where other subtleties/pitfalls 
concerning index D u in the infinite volume setting are also discussed.^ 

To confirm that index D u really is able to capture topological information about 
the continuum gauge field it is desirable to verify that index D u reduces to Q(A) in 
some setting where the former is a well-defined integer from the beginning and the 
problems alluded to above do not arise. To do this, in §5 we consider a setting where 
the continuum field A^x) is in a singular gauge and vanishes outside a bounded region 
of R 4 . We restrict the lattice transcript U of A to a large finite volume hypercube; 

1 Examples of fields with arbitrary non-integer topological charge are easily constructed: View 
R 4 as R 2 x R 2 , let (n, #1) and (r 2 , 82) be polar coordinates on the two copies of R 2 (0j € [0, 2ir]) 
and choose a smooth function A(r) on R + which vanishes in a neighbourhood of and is equal 
to a constant c for r > ro (for some finite ro). Then for any generator T of SU(N) the field 
A = \(n)Td9i ± \(r 2 )Td0 2 has topological charge Q(A) = =Ftr(T 2 )c 2 . This field and its partial 
derivatives are bounded on R 4 . Note that the non-integrality of Q(A) does not contradict the results 

of jl^] since the Euclidean action of the field is infinite. 

2 Even if one allows the values ±00, Q(A) is still not well-defined in general since, e.g., the 



integrand in (1.1) can be oscillatory at infinity. 

3 In |l3| it was erroneously claimed that index I? 17 — * Q(A) in the classical continuum limit if only 

Afj,(x) and the derivatives duA^x) are bounded. We emphasize that this is not true in general: these 

conditions are not enough to guarantee that Q(A) is integer or even well-defined, or that index D u 

is well-defined, cf. the example above. 
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then, with (anti)periodic boundary conditions on the lattice spinor fields, index D u 
is a well-defined finite integer. We show that, in the continuum limit defined by first 
taking the infinite volume limit and then the a — > limit, index D u reduces to Q(A). 
Since index D u and Q(A) are both integers in this setting it follows that they are also 
equal close to this continuum limit, i.e. for all sufficiently large volumes and small 
lattice spacings. This analytic result complements the previous numerical results of 

0- 

The classical continuum limits of the axial anomaly and topological charge in the 
finite volume 4-torus case are not considered here; this case is treated in a separate 
paper |T5|| . 

The organisation of this paper is as follows. §2 provides the relevant background 
material. The material in §3-5 has been described above. (§3: classical continuum 
limit of the index density/axial anomaly; §4: aspects of index/} 17 and its continuum 
limit in the infinite volume lattice setting; §5: classical continuum limit of index D u 
starting in a finite volume lattice setting and taking an infinite volume limit.) The 
conclusions of the paper are summarised in §6. 



2 Background 

In continuum QCD the axial anomaly arises from the triangle diagram in perturbation 
theory ||, and can also be understood non-perturbatively as arising from the jacobian 
for chiral transformations of the fermion fields, regularised by Fujikawa's method 
[|16|1 . In traditional lattice fermion formulations, such as Wilson-Dirac and staggered 
fermions, the axial anomaly arises in a different way: The chiral jacobian is trivial, 
and the anomaly arises instead as a result of the breaking of chiral symmetry by the 



lattice fermion action; see, e.g., [|17], |l8j. (The lattice fermion action needs to break 



chiral symmetry in order to avoid species doubling |L9| 



However, a new perspective is possible for the new lattice fermion actions where 
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the lattice Dirac operator satisfies^ 

-D75 + 75-D = aDj 5 D (Ginsparg- Wilson relation [5]) (2.1) 

D* = 75-D75 (75-hermiticity) (2.2) 

(a=lattice spacing). Besides the overlap Dirac operator (( |2.9| ) below), another so- 
lution is the Dirac operator obtained in the perfect action approach, which is given 
less explicitly via recursion relations The fermion action S = ^J^x^i^Dipix) 



is invariant under a new lattice-deformed chiral transformation of the fermion fields 

1,0 

Sip = %ip , Sip = ^75 (2.3) 

where 

75 = 75(1-^), (2.4) 



It follows from ( |2.1| )- (|2T2| ) that 75 = 1 and 75 = 75 , so 75 determines an orthogonal 



decomposition of lattice spinor fields into "chiral" subspaces just like 75. The fermion 



measure is not invariant under the new chiral transformation (|2.3|) , and the axial 



anomaly can be determined from the corresponding jacobian to be ||[] 

A(x) = ti(%(x,x)) = —iatr("f 5 D(x,x)) (2.5) 

Thus the lattice regularisation with the new lattice fermion actions is providing a 
lattice version of Fujikawa's non-perturbative perspective on the axial anomaly. 

A key question now is whether the axial anomaly fl2.5p for the new lattice fermion 
actions reduces to the usual axial anomaly A con t(x) = 2iq A (x) in the classical contin- 
uum limit. There have been a number of attempts to show this in the perturbative 

4 This and the following structure is implicit in the vector overlap formulation; see, e.g., |po| . 
5 More generally the action is invariant under Sip = 75(1 — (1 — t)aD)ip , S^> — "075(1 — taD) 

Vf G [0, 1]. These transformations all lead to the same axial anomaly ( |2.5| ). 

6 This is the same as the axial anomaly for the vector overlap, which was considered in a special 

case in §10 of B. 
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framework for general lattice Dirac operator satisfying eq.'s (|2.1| )- (|2~2j ) [[5], 23]. A 
very general perturbative argument which encompasses these lattice Dirac operators 
has been given in [p4j] , with the correct continuum anomaly being reproduced pro- 
vided certain general conditions are satisfied. In the specific case of the overlap Dirac 
operator an explicit perturbative calculation was carried out in Subsequently, 



other calculations of the continuum limit of the axial anomaly for the overlap Dirac 
operator were carried out in [10 , vl1,p6||,p7]. These did not involve an explicit per- 



turbative expansion in the gauge field; however they did rely on the gauge field being 
sufficiently small (and hence topologically trivial) so that power series expansions 
of certain operators involving the gauge field could be carried out. The arguments 
also required the gauge field to have a well-defined Fourier transformation, which 
is not always true for general gauge field on R 4 . Thus the problem of providing a 
non-perturbative derivation, valid for general classes of gauge fields (including topo- 
logically non-trivial fields) still remained. To provide such a derivation is one of the 
main purposes of the present paper. Our explicit, detailed derivation of the contin- 
uum limit of the axial anomaly (in §3) complements an abstract, implicit argument in 
the case of general gauge field outlined by M. Liischer in §5 of PBfl . (We should also 
mention that the continuum limit of the axial anomaly in the vector version of the 



overlap formulation was calculated in [29], and in the domain wall formulation in [30], 
although these calculations were based on certain approximations (e.g. linearisation 
of the overlap) and/or assumptions (e.g. weak field, slowly varying field).) 

Another appealing feature of lattice Dirac operators satisfying the GW relation 
( |2.1| ) is that the zero-modes have definite chirality, i.e. the nullspace ker D is invariant 
under 7s : if Dip = then D{^ip) = (-75-D + aD^ 5 D)ip = ||. Thus, at least 
in finite volume settings, index D = Tt^I^d) is well-defined. This is in contrast 
to traditional operators like the Wilson-Dirac operator for which only an approxi- 
mate (non-integer) index can be defined [17]. Furthermore, there is a "lattice index 



theorem" 22 



index D = ~Tr(j 5 D) . (2.6) 



Thus index D = a 4 J2xQ( x ) with the index density 

q{x) = ~tr{ l5 D{x,x)) (2.7) 

related to the axial anomaly fl2.5|) by 

A(x) = 2iq(x) (2.8) 

just as in the continuum. The index of D provides a fermionic definition of the 
topological charge of the background lattice gauge field. In the case of the overlap 
Dirac operator^] 

D = Ul + 1 ,e{H)) , e(H) = ^L = -75 (2.9) 



where 



I# = 75 (D w -r^) (2.10) 

with D w being the Wilson-Dirac operator and r > the Wilson parameter, the index 
formula (|2.6|) givesQ 

index D = --Tr(-JL=) (2.11) 



2 ^ y/m- 

i.e. the spectral asymmetry of — H = —H mo . Since H m has symmetric spectrum and 
no zero- modes for m < |1[ ^ [14] , ( |2.11|) equals minus the spectral flow of H m coming 
from the eigenvalue crossings that occur in the interval < m < mo. This is precisely 
the definition of lattice topological charge arising in the overlap formulation 0. 
(The spectral flow definition of lattice topological charge had earlier been studied 
numerically in This lattice topological charge is well-defined for all lattice 

gauge fields except those for which H mo has zero-modes. After excising this measure- 
zero subset, the space of lattice gauge fields splits into topological sectors labelled 



7 In fact all solutions of (2.1)-( |2^ ) have this "overlap" form: D is a solution if and only if 

D = -(1 + 7 5 e) for hermitian e (= —75) with e 2 = 1 , cf. Neuberger's first paper in Q|. 

8 This is assuming a finite lattice and/or suitable periodicity conditions so that the space of lattice 

spinor fields is finite-dimensional and hence Tr(7s) — 0. The role of the relation Tr(7s) — in this 

setting has been discussed in detail in Q ; related issues are discussed in [[32], [33| . 
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by index/). (This is reminiscent of the way that Liischer's geometrical definition of 
lattice topological charge gives a decomposition into topological sectors after excising 



a measure-zero subset fl35|.) The physical region for the parameter mo, i.e. the values 
for which D is not aflicted with species doubling in the fee fermion case, is < itlq < 2 
JD, H, |{|. As we will see, it is precisely for m in this region that the axial anomaly 
has the correct classical continuum limit. 



3 Classical continuum limit of the index density /axial anomaly 

In this section we consider the continuum limit of the index density q (x) (or equiv- 
alently, the axial anomaly A(x) = 2iq u (x)) for the overlap Dirac operator D u with 
U being the lattice transcript of a continuum SU(n) gauge field A on R 4 . 

Put a hypercubic lattice on R 4 with lattice spacing a. The space of lattice spinor 
fields (functions on the lattice sites with values in C 4 ® C") has the inner product 

(ij 1 ,ij 2 ) = a^M^*M^ (3-1) 

X 

where a contraction of spinor and colour indices is implied. For a given lattice gauge 
field Un(x) the covariant forward (backward) finite difference operators („V~) 
are given by 

V+ip(x) = U^x^x + ae^-tPix) (3.2) 
= i>i x ) ~ U fi( x - ae^tpix - ae^) (3.3) 

where e M denotes the unit vector in the positive /x-direction. Set V M = |(V+ + V~) ; 
this operator is anti-hermitian with respect to the inner product ( |3.1| ) since (V^)* = 
VJ. The Wilson-Dirac operator is now given by 

D w = -^+ r -a(\A) (3.4) 

where ^ = 7 M V M (the 7 At 's are taken to be hermitian so Jv is anti-hermitian), 
A = E M V M + V+ = E M (V+)*V+ = E M (V M )*V; (hermitian, positive) and r > is 
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the Wilson parameter. We then have the hermitian operator 

F m = 7 5 (a^-rm)=7 5 (y + r(|A-m)) (3.5) 

cf. (|2.10|) . To define the overlap Dirac operator 

D = -(l + lb -^L=) H = H mo (3.6) 

some restrictions must be made on the lattice fields (besides excluding the lattice 
gauge fields for which H has zero-modes). This can be done in one of the following 

ways: 

(I) Require j^H < oo. Then the lattice spinor fields form a Hilbert space. Since H is 
bounded (triangle inequalities give \\H\\ < 8 + 8r + rmp), D can be defined via the 
spectral theory for bounded operators on Hilbert space (see, e.g., J3BJ). In this setting 
no restriction on the lattice gauge field U^(x) is required (besides the requirement that 
H has no zero-modes). However, the nullspace ker D need not be finite-dimensional 
in general, so index D is not defined in general: we can have index D = oo — oo. 

(II) Require ip(x) and U^(x) to be periodic in each direction with fixed periodicity 
length L. This is equivalent to the finite volume 4-torus setting with topologically 
trivial gauge fields. The space of such lattice spinor fields is finite-dimensional, and 
H leaves this space invariant, so D can be defined via the usual spectral theory and 
has a well-defined finite index. 

(III) More generally, require ip(x) and U^(x) to be related at opposite boundaries of 
a fundamental domain by a gauge transformation. This is equivalent to the general 
finite volume 4-torus setting where the gauge fields may be topologically non-trivial. 
(This was the setting considered in The space of such lattice spinor fields is 



finite-dimensional and D is again well-defined with finite index as in (II). 

(IV) Restrict the lattice to a finite hyper-cubic region and require that if)(x) is periodic 
(or anti-periodic) at the boundary and U^{x) is trivial at the boundary. For such U 
the operator H leaves the (finite-dimensional) space of these spinor fields invariant, 
so D is defined in the usual way and has well-defined finite index. This setting is a 
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variant of (II); the difference is that we do not hold the edge length fixed when taking 
the (classical) continuum limit. Rather, the limit is taken by first taking the volume 
to infinity with fixed spacing a, and then taking a — > (cf. §5). 

In this paper we focus on the settings (I) and (IV). The finite volume torus setting 
(II)— (III) is treated in a separate paper ||15|| . 

In the following the density 0(x,y) of an operator O on the space of lattice 
spinor fields is defined through Oi[)(x) = a^Y^yOi^x^y^iy). Equivalently, 0{x,y) = 
-tt(%, 0%) where S x (y) = S xy . Note that {^f} is an orthonormal basis for the Hilbert 
space of lattice scalar fields in the infinite volume setting (I). 0(x,y) is a linear 
operator on C 4 <8> C n , and is well-defined when O is a bounded operator on the 
Hilbert space of lattice spinor fields with \ \ < oo. The Cauchy-Schwarz inequality 
gives 

\\0(x,y)\\<- A \\0\\. (3.7) 
cr 

In particular, from ( |2.7| ) and ( |2.9| ) we see that the index density is well-defined: 




Our interest is in the classical continuum limit of this quantity, i.e the a — > limit 
when U^{x) is the lattice transcript of a continuum gauge field A^(x) : 

U^(x) = T exp( / aA^x + (1 - r)ae^) drj 

J 

oo - 

= J2 an dT n ---dTxA li (x,T n )---Ap(x,Tx) (3.9) 

n=0 " , 0<Ti<---<-r n <l 

where A^(x, r) = A ll (x + (1 — r)ae M ) and for simplicity the coupling constant has been 
set to unity. The technical setup within which the a — > limits are taken is as follows. 
We assume an infinite collection of hyper-cubic lattices on R 4 has been specified, 
parameterised by the lattice spacing a, with the following properties: (i) For each 
e > there are only finitely many lattices with a > e. (ii) If ax and a,2 are admissible 
lattice spacings and a 2 < ax then the lattice parameterised by a 2 is a subdivision of 
the one parameterised by ax, i.e. the sites of the latter are contained in the set of sites 
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of the former. (E.g. for each a > we have the lattice with sites aZ 4 and can use 
the collection of these lattices parameterised by a = a , a /2, . . . , a /2 p , . . . for some 
ao > 0.) The property (ii) implies that if x G R 4 is a lattice site in some lattice with 
spacing a then it is also a lattice site in all the other lattices with spacing a' < a. In 
the following, in statements concerning a — > limits the variable x always denotes 
such a point in R 4 ; it is fixed in R 4 and does not change as we go from one lattice to 
another. 

To consider the a — > limit we need the overlap Dirac operator D u = i(l+7 5 -^=) 
to be well-defined for small a, i.e. H should not have zero-modes when a is sufficiently 
small. For technical reasons we will furthermore require that H 2 has a strictly positive 
a-independent lower bound: H 2 > b > for sufficiently small a. The existence of 
such a bound follows from the results of |37| (see also f38|). It was shown there that 



when 1 1 1 — U(p) \ | < e for all lattice plaquettes p, where U (p) is the product of the link 
variables U^x) around p, then there is a lower bound H 2 > b, depending only on e 
and mo, such that for fixed itlq e (0, 2) b > when e is sufficiently small. This result 
generalises to the case of arbitrary mo ^ {0, 2, 4, 6, 8} P5[ . In the present case, where 
the link variable is given by (|3.9|), the plaquette variable has the standard expansion 



1 - U(p x ,w) = a 2 F^ u (x) + 0{a 3 )(x) (3.10) 

and hence 

\\1-U(p)\\ ~ 0{a 2 ) (3.11) 

Strictly speaking this bound requires certain conditions on A^(x). We will discuss 
these below, but for the moment we proceed under the assumption that ( |3.11| ) is valid, 
i.e. that there exists a finite K independent of p and small a such that ||1 — Z7(p)|| < 
a 2 K. Then, for arbitrary e > we have ||1 — Z7(p)|| < e for all sufficiently small a, 
and hence an a-independent lower bound H 2 > b > is guaranteed. 

To evaluate the continuum limit we use an integral representation for l/V H 2 to 
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expand it in a power series as follows.| We first decompose 

H 2 = L - V 



(3.12) 



where 



with 



L = -V 2 + r 2 (±A-m) 2 (3.13) 
V = = r l - 1 ,V, + \[ 1 , )lv }V, u (3.14) 



V, = kvJ + V;),£(V--V+)] (3.15) 

v,u = J[(vj + v;) , (v+ + v;)] (3.16) 

As pointed out in [37], the norms of the commutators of the V^'s are bounded by 
max p \\ \ — U{p)\\. This has the following consequences: First, by ( |3.11| ) 

1 1^| | ~ 0(a 2 ). (3.17) 

Furthermore, 

V(x, y) = in regions where the gauge field is pure gauge. (3.18) 

The former implies that for small a we have < 6/2 where b > is the above- 
mentioned lower bound on H 2 . This in turn implies the lower bound L > 6/2 > 
for the positive operator L in ( |3.13| ). Thus for sufficiently small a the operator L is 
invertible, ||V|| < 1, and we can make the expansion 

H f™ da 1 

H 



H 2 J~oo vr H 2 + a 2 

r°° da / 1 \ / 1 

J-oo ~ { 1 - (L + a 2 )~ 1 V ' \ L + a 



H 



— J2(G a V) k G a . (3.19) 

-°° n k=0 



9 This expansion was used in Kerler's paper [[13] where it was presented as new. In fact it had 
already been given by the present author in fllj] . 
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where G a := [L + u 2 )^ 1 . Note that the 7-matrices in (|3.12 ) are all contained in V. 
Since the trace of 75 times a product of less than 4 7-matrices vanishes, the k = 
and k = 1 terms in ( p,19|) give vanishing contribution to q u (x). On the other hand, 
the terms with k > 3 satisfy the following bound: 

00 71 fc=3 



< 



1 



\H\ 



< a 2 K 3 \\H 



30 da 

00 IX 



fc=3 

00 da 



ife+ii 



^1 



(3.20) 



fc=0 



, 7T (6/2 + a 2 ) 4 _ 

where we have used ( |3.7| ), (|3.11|) and the bounds G a < (6/2 + a 2 ) -1 < 2/6. This is 
0(a 2 ) since the integral and sum are finite and remain so in the a — > limit. Hence 
only the k = 2 term in (|3.19| ) contributes in the a — > limit of the index density 
i.e. 



q U {x) = q U 2 {x) + 0{a 2 ){x) 



where 



qV{x) 



1 r°° drr 

-- — tr[HG a VG a VG a } {x, x) . 

2 J-00 7T 



(3.21) 



(3.22) 



For lattice operators O which are polynomials in we denote by the operator 
obtained by setting U — 1 in ( |3.2j )-(|3~3]). Expanding U^(x) in powers of a via ( |3.9| ) 
gives \ \H — ~ 0(a) and ||L — ~ 0(a). (The rigorous justification of this 

again requires certain conditions on A^(x) , to be discussed below.) The latter implies 
\\G a — G^\ \ ~ 0(a). This can be seen in various ways, e.g. as in ]10|, (v4)], or more 
simply by noting that G a - = G^(L^ - L)G a ; this gives \ \G a - G®\\ ~ 0(a) 
due to the above- noted upper bound G a , G^ < 2/6 which holds for sufficiently small 
a. This allows us to replace H and G a by and G^ in ( |3.22| ) at the expense 
of an 0(a)(x) term. Furthermore, since ||V|| ~ 0(a 2 ) we have V}\\ ~ 0(a 3 ), 

which leads to 1 1 [G^ , V] \ | ~ 0(a 3 ) as follows: The bound ||V^|| < 2 and triangle 
inequalities lead to an a-independent upper bound L < c which allows to expand 

In/ 1 



G a 



c-L 
c+cr 



(3.23) 
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Now, since 



(c- L (0) ) m ,V}\\ < m\\[L {0) ,V}\\ ■ \\c-L\\ m ~ l < m{a 3 K)(c - b/2) m - 1 (3.24) 



we get 



\\[G?\V]\\ < ^E^im+V ' 1,2 



m=0 



(3.25) 



and this is ~ 0(a 3 ) since the sum converges (since < 6/2 < c). Taking this into 
account in ( |3.22| ), it follows from ( p.21[ ) that 

1 r°° da 



q{x) 



2 J-oo 7T 
1 

-tr 



H (0) V 2 



ti[H^V 2 (G^) 3 ](x,x) + 0{a){x 
00 da 1 



2\3 



x, x) + 0(a)(x) 



-^-tr[H {0) V 2 (L {0) )- 5/2 }(x,x) + 0(a)(x) 
16 



(3.26) 



Evaluating the trace over spinor indices we find 
— 3r 



q u (x 



16 



-fiu pa 



tr 



+0(a){x) 



Vf{V u V pa + V V pV a ) + (|A(°) - m )VM(L ( °y)- 5 / 2 \(x,x) 

(3.27) 



where and V^ v are given by ( |3.15| )-( |3.16| ). Noting that [|37| 

[V+ V+]^(x) = (1 - Uip^^U^Uuix + ae^(x + ae^ + ae u ) (3.28) 



and similar formulae for the other commutators, calculations with ( |3.9| ) give 

[Vj,V±]V(aO = (a 2 F M ^(x) + 0(a 3 )(x))^(x ± ae M ± ae y ) 
[Vj,Vj]^(x) = (a 2 i^,(x) + 0(a 3 )(x))ij(x±a eil Tae u ) 



(3.29) 
(3.30) 



These determine the relevant contributions of and V pu in ( |3.27| ). 

We now exploit the fact that the delta-function 8 X on the lattice sites has the 
Fourier expansion in plane wave fields: 

r* d 4 k 



-Tr (2tt) 4 
14 



e" ifc:c/a fc 



(3.31) 



where 4>k(y) '■= e %ky l a . For a general operator O this leads to 

U{X,X) a ^ a 2' U a 2) „.. J , , 



d A k 



- ikx l a {0<p k ){x) 



(3.32) 



In the case where 



O = e 



Vf{V v V pr , + V up V a ) + (1A<°> - m )V, v V pa ){L^)-^ (3.33) 



a calculation using {WM- ^M with (^29|)-(^30l) gives£] 

(0<j> k )(x) = 32n 2 a 4 \(k-r,m )(q A (x) + O(a)(x))Mx) 
where q (x) is the continuum index density (|1.1|), and 

n^cos^^m + E^l-cosA;^) - £ M ^f) 



(3.34) 



A(fc; r, m) 



E« sin 2 fc^ + r 2 (-m + E (1 - cos k p )) 2 



5/2 



(3.35) 



It follows from flOTp and (|3~3p that 

g(x) = I(r,m )q A (x) + 0(a)(x) 



where 



J(r, ml 



-3r 



57T* J-7T 



d k X(k] r, m) 



(3.36) 



(3.37) 



The integral J(r, m) is similar to the integral (A. 17) of | lq| , although the exponents 
and numerical factor are different and the parameter m did not appear there. To 
evaluate it we exploit the symmetries of the integrand (as in JT3) and change variables 
to s v = sin k v to write 

4 



I(r,m)= ( II sign(e /i )J/(r,m,e) 

e M =±l (1=1 



(3.38) 



where 



J(r, m, e) 



-3r r 1 



5^" J-l 



-™+e m (i-wi-^)-£ 



5/2 



(3.39) 



10 Essentially the same calculations were presented as new in [T^]. In fact they had already been 
done in @ (v4)] and Q. 
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This diverges for m = m e = X)u(l — e A*) e {0,2,4,6,8} but is finite for all other 
values of m. It is a constant function of to with a jump at m e ; to see this set 
A(s,m) = -m + Eu(l - W 1 - s l ) , then 



■8tt 2 d r . 



d 4 s 



(3.40) 



/_! ( s 2 + r 2 A 2)5/2 J_ x ^2 +r 2 A 2)7/2 

Inspired by the identity eq. (A. 19) of |18| (which was originally due to Karsten and 
Smit) we rewrite the second integral as 

(l-±E^£)(s 2 + rW) 

-l 

^-"" 2 (3.41) 



(s 2 + r 2 A 2 ) 7 /2 
5r/y s (l + I£ S ^)( S 2 + r 2 A 2 



Integration by parts now gives the first integral in ( 3.40|) except with the opposite 
sign, so ( |3.40| ) vanishes. J(r, to) can now be determined by evaluating J(r, to, e) in 
the limits m — > m e from above and below. For m > m e , after setting m = m e + to 



and changing variables to s v = s u /m in (|3.39| ), we get 



I(r, to, e) 



-3r r 1 /' 



,4 4 

a sm 



—to + 0(m 2 



m— >0+ 



57T J-l/m 
-3 roo 



d 4 ? 



to^s 2 + r 2 (— 1 
-r 



0(to)) 2 ] 5 / 2 
1/2 (3.42) 



J^g2 _|_ r 2^5/2 

For m < m e a similar calculation gives /(r, m, e) = — 1/2. Thus 7(r, to) is independent 
of r > and can now be calculated from ( [3.38 ), leading to the value /(to) in table 
below:Q 



/(to) 



< to < 2 


2 < to < 4 


4 < to < 6 


6 < to < 8 


to g [0, 8] 


1 


-3 


3 


-1 






(3.43) 



By (|3.36 ) we now have 



q u (x) = I(m )q A (x) + 0(a)(x) 



(3.44) 

11 The evaluation of this integral was given in the first version M, (vl)] of this paper, and also (in 
more detail) in |27||. 
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Thus for < m < 2 the index density q (x) does indeed have the correct classical 
continuum limit in the infinite volume setting. 

The 0(a p ) bounds used in the calculations above can be established by standard 
calculations when the continuum field A^x) is smooth and ||^4 M (a;)|| , and 
| \d ll d v A p {x) 1 1 are all bounded on R 4 (the details can be found in the appendix of 
). However, the result q(x) = q (x) + 0(a)(x) (and its generalisation fl3.44|) when 



mo (jL [0, 2]) holds in more general cases as well: It suffices that the plaquette variable 
of the lattice transcript of A satisfies 1 1 1 — U(p) 1 1 < e when the lattice spacing a is 
sufficiently small, with e small enough that a strictly positive lower bound H 2 > b > 
is guaranteed (as discussed earlier). Then D u is local in the gauge field [j37| , leading 
to 

\q u (x)-q"(x)\ ~ 0(±e-^ a ) (3.45) 

where U is the lattice transcript of a smooth continuum field A which coincides with 
A in a neighbourhood of x and vanishes outside a bounded region of R 4 . Since such 
A and its partial derivatives are automatically bounded on R 4 , we have q u (x) = 
q A (x) + 0(a)(x), and this together with ( |3.45|) gives q(x) = q A (x) + 0(a)(x). 

A more detailed justification of the preceding is as follows. Let A be an arbitrary 
smooth continuum gauge field on R 4 . Pick a smooth function X(y) on R 4 which is 
equal to 1 in the box ELKu — de^, x^ + de M ] C R 4 (d > 0) and which vanishes outside 
a bounded region of R 4 . To prove the claims above, it suffices to establish (|3.45 ) 



for the case where A^y) = X(y)A fl (y). For this we exploit the fact [[57|] that there 
is a power series expansion l/\/lP = KJ2kLot k Pk(H 2 ) where Pfc(-) is a Legendre 
polynomial of order k; \\Pk{H 2 )\\ < 1; t — e~ e ; the constants k, 9 > depend only 
on the (strictly positive) lower and upper bounds on H 2 P7[| . Similarly 1/V H 2 = 
K J2T=o t k Pk{H 2 ) where H is obtained from H by replacing the lattice gauge field 
U (the lattice transcript of A) by the lattice transcript U of A. (This requires H 2 
and H 2 to have a lower bound b > 0. This is guaranteed in both cases when a is 
sufficiently small: In the first case we have ||1 — Z7(p)|| < e by assumption, while 
in the second case ||1 — ~ 0(a 2 ) since A^(x) and its partial derivatives are 
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automatically bounded on R 4 .) Since H only couples nearest neighbour sites, P k (H 2 ) 
can only couple the site x to itself via a site outside of rLKu — de^, x p + de p ] C R 4 if 
k > 2(d/2a). Therefore [P k (H 2 )}(x, x) = [P k (H 2 )}(x,x) when k < d/a, and it follows 
that 

( x ,x)-- m ,{x,x)\\ < k t k \\[P k (H 2 )](x,x)-[P k (H 2 )](x,x)\\ 

k>d/a 



oo -r 



K _ e )- A e-°^ (3.46) 
1 — e u ' a 4 



This together with the fact that H is ultra- local gives Q3.45Q . 



The preceding observations allow us to conclude that q u (x) = q A (x) +0(a)(x) in 
some cases where A^x) diverges at infinity, or has singularities. An example of the 
former is a topologically non-trivial gauge field on the 4-torus: these can be viewed 
as gauge fields on R 4 satisfying a periodicity condition 

A p (x + Le u ) = n(x, v)Ap(x)Sl{x, z/)" 1 + Q(x, u)d,Sl(x, v)' 1 (3.47) 

where Q(x,v) , v — 1,2,3,4, are the S l f/(n)-valued monodromy fields which specify 
the principal SU(n) bundle over T 4 . Fields A^x) satisfying ( |3.47| ) diverge at infinity 
in the topologically non-trivial case. However, the requirement ||1 — C/(p)|| < e for 
small a still holds in this case since (i) the lattice transcript of A^x) satisfies 

Unix + Le u ) = Q(x, v)U^(x)VL(x + ae u , z/)" 1 (3.48) 

and (ii) ||t^t(|/)|| = 1 for all y y \i since U^y) is unitary. These imply that if x' p = 
x p + Ln p e p then ||1 — U{p x ^ u )\\ = ||1 — U{p x > >(il> )\\ for for all n G Z 4 . It follows 
that ||1 — Z7(p)|| ~ 0(a 2 ) since this is true for plaquettes in the fundamental domain 
[0, L] A C R 4 because A p (x) and its partial derivatives are automatically bounded in 

An example where the requirement 1 1 1 — U(p) 1 1 < e for small a is satisfied when 
Ap(x) has a singularity is the following: Consider a gauge field in a singular gauge, 
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and such that A^(x) is pure gauge in a neighbourhood of the singularity and vanishes 
outside a bounded region of R 4 . (Examples of topologically non-trivial gauge fields 
of this type are readily constructed, cf. §5.) We choose the lattices on R 4 such that 
the singularity of A^ix) never lies on a link of any lattice (cf. §5). Then the lattice 
transcript U^x) is well-defined for all lattices. For lattice plaquettes contained in the 
neighbourhood of the singularity we have U(p) = 1 since A^(x) is pure gauge. On the 
other hand, A^x) and its partial derivatives are bounded outside the neighbourhood 
of the singularity (since the field vanishes outside a bounded region). Hence ||1 — 



4 Aspects of index D and its continuum limit in the infinite 
volume setting 

In contrast to q u (x), the index of D u is a problematic quantity in the infinite volume 
setting. It is not well-defined a priori by 



since the null-space ker D may be infinite-dimensional. The same is true for the index 
of the continuum Dirac operator^ for general gauge field A^{x) on R 4 . In the latter 
case the index exists and is equal (by the index theorem) to Q(A) = J R4 d A x q (x) 
provided A^(x) converges sufficiently quickly to a pure gauge configuration at infinity. 
A natural conjecture in the lattice setting is therefore that index D u , given by ( |4.i|) , 
exists and is equal to a 4 J2 X Q U i x ) when the lattice gauge field U^x) converges suffi- 
ciently quickly to a pure gauge configuration at infinity. To prove this is a challenging 
problem though, and in this section we will only make some first steps towards it. We 
investigate what happens if q u (x) is replaced by a 4 J2 X Q U i x ) m the continuum limit 
calculation of the preceding section; this will indicate how a situation where £7 M (x) is 
pure gauge at infinity can result in the infinite sum J2 X Q U ( x ) being convergent and 
° 4 J2 X q u (x) — » Q(A) in the classical continuum limit. However we do not show that 



U( P )\\ - 0(a 2 ). 




(4.1) 
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index D u exists and is equal to a 4 J2x Q ( x ) i n this case; this remains as a problem 
for future work. 



Using the expansion (|3.19|) we get 



| i-OO (iff 00 

a 4 E^C*) = --a 4 £ \ — Y.^\m«V) k G«\M (4.2) 

x 1 x J -°° n k=0 

As before, the k = and k = 1 terms give vanishing contribution. Obviously, to get 
° 4 J2 X q U (x) f-Rid^xq (x) = Q(A) we will need the part with k > 3 to vanish in 
the a — > limit as in §3. This cannot be expected to happen in general since the 
infinite sum J2 X might not even be convergent for k > 3 part of (|4.2| ). However, as 
we will now show, the sum over lattice sites does converge when U^{x) is pure gauge 
outside a bounded region 1Z of R 4 . Furthermore, we will see that when U^(x) is 
the lattice transcript of a smooth continuum field A^(x) which is pure gauge outside 
of 1Z the k > 3 part in ( f4.2|) does indeed vanish in the a — > limit and we get 
a 4 J2x q u {x) —* Q{A). We can assume that 1Z is a 4-dimensional box [-L, L] 4 . Let TZ^ 
denote the larger box [— L — d, L + d} 4 (d > 0). The strategy is to split up the sum over 
lattice sites in (|4.2|) into a finite sum Y^ x eR d (' ' ') an d an infinite sum X^eR 4 -7£ d (' ' ') 
where the sums are over the lattice sites in IZd and R 4 — TZd respectively. The first 
sum is easy to deal with after noting that a 4 Z^e^G ' ') * In d ^ x (' ' ') f° r a ~^ ® 
and using the fact that the volume V(7Zd) is finite. We will deal with the second sum 
by exploiting the facts that G a (z, x) is (exponentially) local and V(y, z) is ultra-local 
and vanishes in regions where U is pure gauge. 

By the same calculations as in ( |3.20|) we see that the k > p part of the summand 
in 04.2|) has a bound 

k=p 



/OO rjfj °° 
— J2^[H(GaV) k G a ]{x,x) < a 2 ^K p (4.3) 
-00 71 , 



(p > 2) for some constant K p independent of a and x. Choose an e > 0, then for 
sufficiently small a we have a 4 J2x&Ti d 1 — VCR-d) + e = V(7Zd) t where V(7Zd) = 
(2(L + d)) 4 is the volume of Rd- It follows that, for small a, 

da 00 
xeiZd k=p 
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a 4 



/OO riff 00 
— J2tv[H(G a V) k G a ](x,x) < a 2 ^K p V{K d ) t (4.4) 
-co 71 , 



This vanishes in the a — > limit for p = 3, and does not diverge for p = 2. Hence 
when the sum over x in (|4.2| ) is restricted to the lattice sites in IZd the k > 3 part 
vanishes in the a — > limit as required, and the k = 2 part remains finite. We now 
consider the summand in (|4.2|) for x outside of TZd- We use 

(4.5) 



[H{G a V) k G a ](x, x) = a 8 Y l [H(G tr V) k - 1 G cr ] (x, y)V(y, z)G a (z, x) 



l).Z 



By (|3.18|) V(y, z) = for y, z ^ 1Z so the sum over y and z in ( |4.5| ) can be restricted 
to the lattice sites in 1Z. We now apply a version of the locality argument of to 
Ga : For small a we have a-independent bounds < b/2 < L < c and get a power 
series expansion of G a as in |3.23| , 



G n 



c + ° r2 ^o Vc + or2 



(4.6) 



which converges since ||^rf || < ^zr0 < < ~l~ = t < 1. Since L only couples nearest 
neighbour and next to nearest neighbour sites, L k (x, y) = when J2u \ x ^~Un\/ a > 2fc. 
It follows that 



I G & (x, 



< 



1 



c + a 2 



E 



C — L \ k 



c + a' 



< \\GJ\t 



oo 

J2t k /a- 

fc=0 



< \\G a \\ — exp(--J2\Xf,-yf,\/a 
a z n 



(4.7) 



where t = e 8 = ° b ' 2 and k = 1/(1 — t). Applying this to G a (z, x) in (|4.5| ) gives 
\\G„{z,x)\\ < ||G CT ||4exp(-|5:(K|-L)/a) (4i 



since |x M — 2^1 > |x„| — L when z ElZ and x is outside of 7^. This leads to a bound 



on the part of (|4.2|) where the sum is restricted to the lattice sites in R 4 — 1Z 



E 



roo r] n 00 

_ Y.^\H(G a V) k G a ]{x,x) 

xeR, 4 -1Z d J 00 71 k=0 

da 



/OO fjfj 
— £ ti[H{G (r V) k - 1 G a ] (x, y)V{y, z)G a {z, x) 

< « 4 E E/ — E^ll^vf-^iiiiFiiHG^^ii 

xeR 4 -7e d y,ze7?- 00 71 k=i 
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< 



1 r°° rlrr 00 

' —Y,M\H(G*V) k G a \ 

00 71 k=2 



,1 




X 



9 

x a . 

^ n 4 V 9 
xeR 4 -K d " 



' E 4«p(-?HW-i)/«) (4-9) 



Since ||V|| ~ 0(a 2 ) the integral is ~ 0(a 4 ). For small a we have a 8 J2 y ,zeii^ < 
(V(7?.) e ) 2 . Finally, the sum is bounded by 



n 



/■oo / H \ / 7 \ 4 



Thus (|4.9|) vanishes as 0(^e~ p / a ) for a — >• 0. This completes the demonstration 
that the sum over lattice sites in ( |4.2j ) is convergent when U is pure gauge outside a 
bounded region, and that the k > 3 part vanishes in the a —>■ limit when [/ is the 
lattice transcript of a continuum field A which is pure gauge outside a bounded region. 
Furthermore, the preceding shows that in this case the non-vanishing contribution to 
a 4 J2 x Q U ( x ) i n the a — >• limit comes from a 4 J2 x &n d I2 ( x ) where the subscript "2" 
refers to the k = 2 part of (|4.2|) . It is a straightforward technical exercise to show that 
this reduces to j n d A xq A {x) in the a — >• limit along the same lines as the argument 
for q u (x) — > in §3 (we omit the details). By Stokes theorem this reduces to 

the integral of the Chern-Simons term over the boundary dTZj, which calculates the 
winding number of the map from dlZd = S 3 to SU(n) corresponding to the pure 
gauge configuration, and this number is precisely Q(A). 

Clearly the preceding arguments can be generalised from the case where U and 
A are pure gauge outside a bounded region of R 4 to cases where the fields converge 
sufficiently quickly to pure gauge configurations at infinity. To determine a precise 
criterion for what is "sufficiently quickly" is a non-trivial problem though, and is 
beyond the scope of this paper. In the continuum the criterion is that F^x) should 



vanish quickly enough at infinity to be square- integrable |12| . We speculate that a 
similar condition will suffice in the lattice setting. 
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5 Classical continuum limit of index D u starting in a finite 
volume lattice setting and taking an infinite volume limit 



In this section we consider index D u in the lattice setting (IV) of §3. The volume 
is finite, hence the space of lattice spinor fields is finite-dimensional and index D u 
is well-defined from the beginning. The boundary condition on the lattice gauge 
field is that it is trivial at the boundary (i.e. the Dirichlet condition). Then the 
covariant finite difference operators map both the spaces of periodic and anti- 
periodic lattice spinor fields to themselves, and in both cases satisfy ( V^)* = — , 
so H is hermitian and D u is well-defined. We can therefore take either periodic or anti- 
periodic boundary conditions on the lattice spinor fields. We consider the case where 
U is the lattice transcript of a continuum gauge field A^(x) in a singular gauge and 
vanishing outside a bounded region of R 4 . The singularity allows A to be topologically 
non-trivial (i.e. Q(A) 7^ 0) while the latter requirement ensures that U^x) is trivial 
outside of a bounded region of R 4 and therefore satisfies the Dirichlet boundary 
condition when the volume is sufficiently large. We will show that index D u reduces 
to the continuum topological charge Q(A) in the classical continuum limit defined by 
first taking the infinite volume limit and then taking a — > 0. Since index D u and Q(A) 
are both integers, it then follows that index D u = Q(A) for all sufficiently large lattices 
with sufficiently small lattice spacings. This demonstrates analytically that index D u 
is able to capture the topological data of the continuum field A (complementing 
previous numerical results which we discuss later). The present setting is arguably 
the simplest in which such a result can be analytically demonstrated. The same result 
holds in the finite volume torus setting |n|; however, the argument is less simple than 



what follows.0 

Let A^{x) be an SU(n) gauge field on R 4 which may be singular at the origin but 
is smooth everywhere else and vanishes outside a bounded region. Furthermore we 
require that A is pure gauge in a neighbourhood of the origin. Examples of such fields 



12 In fact the result derived in this section can also be derived quite straightforwardly as an appli- 
cation of the result for the finite volume torus case |^(| . 
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are readily obtained as follows. Take a smooth map <fi : S 3 — > SU(n) with degree 
Q and define <f> : R 4 — {0} — >■ SU(n) by 4>{y,t) = 4>{y) , y E S 3 , t <E R + , where we 
are identifying R 4 — {0} with S 3 x R + in the obvious way. Choose a smooth real 
function X(x) on R 4 equal to 1 in a neighbourhood of the origin and vanishing outside 
a bounded region. Then the field 

A^x) = Xix^d^ix)- 1 (5.1) 

has topological charge Q, is singular at the origin, pure gauge in a neighbourhood of 
the origin and vanishes outside the bounded region. 

We take the finite hyper-cubic lattice in R 4 to have spacing a and sites 

{ x = a (m - 1/2 , n 2 - 1/2 , n 3 - 1/2 , n 4 - 1/2)} \ n E Z 4 - N < < N } (5.2) 

Note that the origin G R 4 never lies on a lattice link, so the lattice transcript U^x) 
of Afi(x) (given by Q3.9|) ) is well-defined for all a and N. 

Our aim in this section is to show that, with either periodic or anti-periodic 
boundary conditions on the lattice spinor fields, 

lim lim index D u = I(m )Q(A) (5.3) 

a^O N—>oo 

where /(mo) is given by the table ( |3.43| ) (in particular I(m ) = 1 for the physical 
values < m < 2). The techniques used to derive this have already been developed 
in the earlier parts of this paper, so we will be economical with the details in the 
following. 

Proof of j\5.3j ). We start from the formula 

1 H 

index D u = a 4 ^2 q u (x) , q u (x) = — tr (x, x) (5.4) 

x 2 VH 2 ' 

where the sum is over the lattice sites ( |5.2j ). Choose L sufficiently large so that A^(x) 
vanishes for all x outside of 1Z = [-L, L} 4 . Choose d > and set IZj = [—L — d, L + d} 4 
as before. By the argument at the end of §3 we have 1 11 — U(p) \ \ ~ 0(a?) . This leads to 
a lower bound H 2 > b > for small a, with b independent of both a and N. Therefore 
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the previous bounds such as ||V|| < a 2 K and L > b/2 continue to hold in the present 
setting and can be chosen independent of N. We henceforth restrict to the (sufficiently 
small) lattice spacings a for which these bounds are satisfied. Furthermore, for given 
a we restrict to the (sufficiently large) iV's for which aN > L + d so that the lattice 
(specified by ( |5.2p ) covers TZd- The N — > oo limit of the part of (|5.4j) where the sum 
is restricted to sites outside of TZd has the bound 



lim 

N— »oo 



a 4 £ q u (x) < a 4 ]T \q u {x)\ (5.5) 



where the sum on the right-hand side is over all lattice sites outside of R 4 — TZd in 
the extension of the lattice to the whole of R 4 . This vanishes in the a — > limit by 
a similar argument to the one based on ( f4.9|) — (4.10) in the previous section. Thus in 



deriving (|5.3|) it suffices to restrict the sum in (|5.4j) to the lattice sites in TZd- 

Now choose V > small enough that [— L',L'] 4 is contained within the neigh- 
bourhood of the origin in R 4 where A^(x) is pure gauge. Choose e > so that e < U 
(e.g. e = L'/2) and set TZ e = \—V + e, U — e] 4 . By a straightforward adaptation of 
the argument based on ( O )-( CTO ) we find that the part of ( |5.4j ) where the sum is 
restricted to the lattice sites in TZ e vanishes in the a — ► limit after taking N — >• oo. 
Thus in deriving ( |5.3|) it suffices to restrict the sum in ( |5.4j) to the lattice sites which 
are contained in TZd — TZ e , i.e. it suffices to show 

lim lim a 4 Y q u (x) = I(m )Q(A) . (5.6) 

x&TZa — rCe 

Using the locality of q (x) in the gauge field (cf. §3), arguments similar to those in 
§3 lead to an analogue of ( |3.36| ): In the present setting the lattice delta-function has 
the Fourier expansion (analogue of ( |3.31j )) 



-Tr<k u <n \ I 



where (f>k(y) '■= e %ky l a is the lattice plane wave field with momentum k ; the domain of 
the lattice 4- momentum (i.e. the summation domain in (|5.7|) ) in the case of periodic 
or anti-periodic boundary conditions is 

periodic: e ^{-N, -N + 1, . . . , N - 1} (5.8) 
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anti-periodic: G ^{-N + 1/2, -N + 3/2, . . . , N - 1/2} (5.9) 

and A 4 k = 7r 4 /(2iV) 4 = the "volume per k" in lattice momentum space. Using this 
we find, analogously to ( |3.36| ), 

q u (x) = I(r, m , N)q A (x) + 0{a){x) (5.10) 

where 

I(r,m,N) = — — y A 4 k\(k;r,m) . (5.11) 

— n<k u <7r 



with X(k] r, m) given by ( 3-35 ). The iV-dependence is in the summation domain: ( |5.8|) 
or ( |5.9|) . In the N — > oo limit A 4 k — > <i 4 £; and we see from (|5.8|) — (|5.9|) that in both 
the periodic and anti-periodic cases limjv^oc I(r,m,N) = I(r,m) = I(m) (cf. §3 for 
the last equality). It follows that 

lim lim a 4 V q u (x) = I(m )[ d 4 xq A (x). (5.12) 

By Stokes theorem this reduces to the difference of the integrals of the Chern-Simons 
term over the boundaries &Rd and &R e . The former vanishes since dTZd lies in the 
region where A^(x) vanishes. The latter gives the the winding number of the map 
from &R e = S 3 to SU (n) corresponding to the pure gauge configuration, and this is 
precisely Q(A) in the present singular gauge case. This completes the proof of fl5.3Q. 

Comparison with numerical results. Our result ( |5.3|) implies that for < niQ < 2 the 
index of D u reproduces the continuum topological charge Q(A) when N is sufficiently 
large and a is sufficiently small, and that the index reduces to I(mo)Q(A) for general 
mo ^ {0,2,4,6,8}. This is something which can be investigated numerically, and 
in fact numerical investigations in closely related lattice setups have already been 
carried out. In finite volume settings such as the present one (where the space of 
lattice spinor fields is finite-dimensional) the index of D u equals the spectral flow of 
— H m as m increases from any negative value to m Q. This follows from ( |2.11|) : 
index D u = — §Tr(n||) (recall H = H mo ) and the fact that the spectrum of H m is 
symmetric and without zero for m < [[[], 0, [L4[] . In 0| the spectral flow of H m was 
studied numerically for various smooth SU(2) instanton fields on a finite lattice.0 
13 The spectral flow of H m for various lattice gauge fields was first studied numerically in Q. 
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One of the situations considered in [13 was instanton field in a singular gauge with 
anti-periodic boundary conditions on H m . I.e. anti-periodic boundary conditions on 
the lattice spinor fields, and the lattice transcript of the instanton field modified at 
the boundary links to make it anti-periodic, so that H m is a well-defined hermitian 
operator on these spinor fields. This is a minor modification since the localisation 



radius of the singular gauge instanton in [14] was well within the region covered by 
the lattice, i.e. it was almost vanishing at the boundary of the lattice. This is similar 
to the setup that we have considered in this section: the singular gauge instanton 
can be approximated by a gauge field A^(x) with singularity which is pure gauge in a 
neighbourhood of the singularity and vanishes outside a bounded region. The results 
of (|TJJ] are compatible with ours: They numerically determined the crossings of the 
origin by low-lying eigenvalues A(m) of H m for < m < 2 and found precisely one 
crossing, which occurred reasonably close to zero (at m « 0.5), and the slope of A(m) 
at the crossing was such that the spectral flow of — H m from this crossing agreed with 
the topological charge of the instanton (cf. Fig. 7 of fll4|). This can be equivalently 
expressed as index D u = Q(A) for 0.6 < mo < 2. This is an approximative numerical 
confirmation of our analytic result in this section. Numerical studies of the spectral 
flow of H m in other situations in |34] , |i~4f , or equivalently, numerical studies of index D u 
as a function of mo f3^j, are also compatible with the the classical continuum limit 
result for index D u that we have shown for the particular situation considered in this 
section (and which has also been analytically shown in the finite volume torus case 
in |Tf5|): For "sufficiently smooth" lattice gauge fields it was found that (i) crossings 
of the origin by eigenvalues A(m) of — H m occur close to m = 0, 2, 4, 6, 8 ; if the net 
spectral flow from crossings close to m = is Q then the net spectral flows from 
crossings close to m = 2, 4, 6, 8 are, respectively, — AQ , 6Q , — AQ , Q , and (ii) if U 
is the lattice transcript of a continuum field A then Q = Q(A). (A situation where 
this does not hold is when U is the lattice transcript of an instanton field in a regular 
gauge and anti-periodic boundary conditions are imposed on H m ||14|| . This is to be 
expected though, since in this case the instanton field is approximately pure gauge 
at the boundary of the region covered by the lattice, and is therefore not close to 
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satisfying anti-periodic boundary conditions.) 
6 Conclusion 

In this paper we have rigorously verified that the index density/axial anomaly for the 
overlap Dirac operator has the correct classical continuum limit in the infinite volume 
lattice setting, and in a finite volume setting where the the continuum limit involves 
an infinite volume limit. The only condition on the continuum field A^(x) required 
to establish this result is that the plaquette variable of its lattice transcript satisfies 

|[l-Z7(p)|| < e(m ) (6.1) 

for sufficiently small lattice spacing a, with e(mo) small enough to guarantee the 
existence of a lower bound b(mo) > on H 2 |37[ |38|, [39]], which in turn guarantees 
the locality of D u in the gauge field [j3~?| . The locality plays a central role in our 
arguments. The condition ( |6.1|) is automatically satisfied for small a when A^(x) 
and its first and second order partial derivatives are bounded on R 4 , since it can 
then be shown that ||1 — ~ 0(a 2 ). However, using the locality property 

we have seen that the condition is satisfied for more general gauge fields, including 
some cases where the field has a singularity or diverges at infinity. In particular our 
results cover the case of topologically non-trivial gauge fields on R 4 . This case was 
was not covered by previous works on this topic (cf. the discussion in §2). Our 
approach was to decompose the operator H 2 appearing in the overlap Dirac operator 
D = ±(l + 7s -^=) into H 2 = L-V where all the Kamma-matrices are contained in 
V and ||V|| ~ 0(a 2 ). Then, starting with an integral representation for the inverse 
square root, we expand 1/v ' H 2 in powers of V. It was found that only the term of 
order 2 in V contributes to the index density q u (x) = — ^ti^^(x,x) in the classical 
continuum limit, and the contribution from this term was shown to reduce to the 
continuum density q A (x) (for < m < 2). 

The expansion of the overlap Dirac operator in powers of V is of interest in its 
own right and may have other useful applications. Unlike an expansion in powers 
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of the gauge field, such as the one given in [Po| , which requires the gauge field to 
be small (i.e. link variables close to 1), our expansion only requires a bound of the 
form (|6.1|) to hold, i.e. plaquette variables close to 1, and can therefore be used for 
topologically non-trivial gauge configurations. 

Subtleties/pitfalls which arise when considering the index of the overlap Dirac 
operator in the infinite volume setting were discussed (§4). In this case the index 
is generally ill-defined a priori. Formally we have index D u = a 4 J2 X Q U ( X ), an d we 
pointed out that the previous result q u (x) — > q A (x) for the index density does not 
imply a^Y.xQ 11 ^) ~^ fn,* q A (x) = Q(A) in general (as was erroneously claimed 
in jT3|): The argument was seen to break down due to the infinite sum over lattice 
sites in J2xQ U ( x )- We showed however that a 4 J2x 1 U (%) —> J R 4 d A x q A (x) = Q(A) in 
the case where A^(x) is pure gauge outside a bounded region of R 4 . This is a first 
step towards showing the following conjecture: index D u is well-defined and equal to 
a 4 J2 X q u (x) when the lattice field U^{x) converges sufficiently quickly to a pure gauge 
configuration at infinity, and in this case index D u — > Q(A) in the classical continuum 
limit. This conjecture (the continuum analogue of which is known to be true |p"2f ) 
remains as an interesting problem for future work. 

In order to verify that the index of the overlap Dirac operator really is able to cap- 
ture topological information we considered (in §5) a finite volume lattice setting where 
index D u is well-defined from the beginning, namely when U is the lattice transcript of 
a continuum field A^x) which is in a singular gauge and vanishes outside a bounded 
region. (For technical reasons we also required A^x) to be pure gauge in a neighbour- 
hood of the singularity. We described how such gauge fields can be readily constructed 
from maps : S 3 — > SU(n).) We showed that lim a ^ hniAr^oo index D u = Q(A) , 
where N — > oo is the infinite volume limit. This result complements a previous 
numerical result obtained in a similar setting in [I3|. 



For the finite volume 4-torus case the results lim^o q u (x) = q A (x) and lim a ^ index D l 
= Q(A) are established for general gauge fields in a separate paper In that case 
the result for q u (x) immediately implies the result for index D u since the volume is 
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finite. However, the finite volume causes additional techinical complications in the 
calculation of \im a ^ q u (x) ; we refer to |T5| for the details. 

Finally we mention that, while the settings we have considered were all in 4 
dimensions, everything generalises straightforwardly to arbitrary even dimension 2m, 
with gauge group SU (n) for 2m > 4 and gauge group U (1) in the 2-dimensional case. 
In the latter case our result q u (x) — > q A (x) confirms a numerical result in §10 of 
where the axial anomaly A u (x) = 2iq u (x) was numerically determined and compared 
to the continuum axial anomaly for a particular topologically non-trivial U(l) gauge 
field on the 2-torus.0 
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